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Abstract 

We calculate the B ^ D and B —> D* form factors at zero recoil in Staggered Chiral Perturbation 
Theory. We consider heavy-light mesons in which only the light (u, d, or s) quark is staggered; 
current lattice simulations generally use a highly improved action such as the Fermilab or NRQCD 
action for the heavy (b or c) quark. We work to lowest nontrivial order in the heavy quark 
expansion and to one-loop in the chiral expansion. We present results for a partially quenched 
theory with three sea quarks in which there are no mass degeneracies (the “H-U-l” theory) and 
for a partially quenched theory in which the u and d sea quark masses are equal (the “2-|-l” theory). 
We also present results for full (2-1-1) QCD, along with a numerical estimate of the size of staggered 
discretization errors. Finally, we calculate the finite volume corrections to the form factors and 
estimate their numerical size in current lattice simulations. 

PACS numbers: 11.15.Ha, ll.30.Rd, 12.38.Gc 
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I. INTRODUCTION 


The CKM matrix element \Vcb\, which places an important constraint on the apex of the 
CKM unitarity triangle through the ratio |Vub/Kfe I, can be determined from experimental 
measurements of exclusive semileptonic 5-meson decays combined with theoretical input. 
Because experiments measure the product (5(1) ■ iKbl)^, where 5(1) is the B ^ D or 
B ^ D* hadronic form factor at zero recoil, the precision of \Vcb \ is limited by the the¬ 
oretical uncertainty in 5(1). Although, in principle, both form factors can be calculated 
nonperturbatively using lattice QCD, in practice, direct calculations of the B ^ D and 
B ^ D* hadronic matrix elements are plagued by large statistical and systematic errors. 
Hashimoto et al. therefore proposed a method for calculating 5(1) on the lattice using double 


m. 


ratios of matrix elements in which most of the statistical and systematic errors cancel 
This method provides the key theoretical ingredient necessary to allow a precise lattice de¬ 
termination of the 5 —> 5 or 5 —> 5* form factors, and hence a precise determination of 
iKfel- 


li> 


and 


Recently was calculated using 2-1-1 flavors of dynamical staggered quarks 

the analogous calculation of Bb^d* will be done in the near future. Because staggered 
quarks are computationally cheaper than other standard fermion discretizations, staggered 


simulations offer the lightest dynamical quark masses currently available 


0 . 


This result 


therefore has a smaller systematic error associated with chiral extrapolation than previous 


quenched results P|. It is known, however, that the 0{a?) discretization errors associated 
with staggered fermions are numerically significant in current lattice simulations and must 
also be accounted for in the chiral and continuum extrapolation of staggered lattice data . 
This procedure is well-established in the light meson sector: use of staggered chiral per¬ 


turbation theory 


Q,Q,y, s 


functional forms for extrapolation of staggered lattice data 


has allowed precise determinations of light meson masses, meson decay constants, and even 


0 


quark masses p . Staggered chiral perturbation theory was recently extended to heavy-light 


mesons (in which only the light quark is staggered) by Aubin and Bernard HI, and has 
been successfully used in the extrapolation of the 5-meson decay constant (ll|. 

In this paper we use heavy-light staggered chiral perturbation theory to calculate the 
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B ^ D and B D* form factors at zero recoild The resulting functional forms can then 
be used to extrapolate staggered lattice data to the continuum and to the physical pion 
mass. Accounting for staggered discretization errors in this way is essential for a precise 
lattice determination of these form factors, and consequently of \Vcb \- 


This paper is organized as follows. We review staggered chiral perturbation theory for 
heavy-light mesons in Section IHl We then calculate the B ^ D and B —>■ D* form factors 
at zero recoil for a l-l-l-l-l PQ theory, a 2-1-1 PQ theory and full (2-1-1) QCD in Section HTH 
Next, in Section HVl we plot the B ^ D and B ^ D* form factors using reasonable values for 
the quark masses and lattice spacing both with and without the taste-symmetry breaking 
contributions. The dramatic change in the shape of the B D* form factor illustrates 
the necessity of accounting for taste-breaking in the continuum and chiral extrapolation 


f. Q 


of staggered lattice data. In Section |V| we use the method of Ref. M to calculate the 
hnite volume corrections to the B ^ D and B ^ D* form factors. We then estimate 
the numerical size of these hnite volume corrections in current lattice simulations; we hnd 
them to be very small - only one part in 10^. In Section IVTl we conclude. The Appendix 


contains additional 
conventions of Ref. 


brmulae necessary to understand our form factor results. It follows the 


II. STAGGERED xPT WITH HEAVY-LIGHT MESONS 

In this section we review staggered chiral perturbation theory (S^PT) for heavy-light 
mesons, which was developed in Ref. [r| . 

We hrst construct the portion of the heavy meson chiral Lagrangian that only contains 
light quark helds. We consider a partially quenched theory with n havors of staggered light 
quarks. The detailed construction of the leading-order effective staggered chiral Lagrangian 
is given in Ref. [3|; we simply present the results that are necessary for the calculation of 
the B ^ D and B ^ D* form factors. 

We assume that spontaneous symmetry breaking of the S't/(4n) chiral symmetry by the 
vacuum, 

SU{4n)L X SU{4:n)R SU{4:n)v , (1) 


^ Note that a subset of our results was presented in Ref. [l^ . 
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leads to 16n^ — 1 pseudo-Goldstone bosons, which we will generically call pious, that can be 
collected into an SU (4n) matrix: 

S = exp(^<l>//). (2) 


The matrix, <h, which contains the pion helds, is traceless with 4x4 submatrices: 

U TT+ K+ ■■■ 

TT- D ■■■ 

K- S ■■■ 

16 

U = etc. (4) 

where the SU{A) generators, 

Ts — {^5i j (5) 

are Euclidean gamma matrices and is the 4x4 identity matrix. The leading order pion 
decay constant, /, is approximately 131 MeV. Like the pion matrix, the quark mass matrix 
is also 4n x 4n, but it has trivial taste structure: 

0 0 

0 nidi 0 
0 0 nisi ■ ■ ■ 

Under chiral symmetry transformations. 






s ^ 


(7) 

M ^ 

LMR^, 

(8) 

L e S[/(4n)L, 

R e SU{4n)R. 

(9) 


The standard SyPT power-counting scheme is: 

pI/K ~ "^^Mqcd ~ , 


( 10 ) 
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so the lowest-order, 0{p^,mq,a^), staggered chiral Lagrangian is^ 


rLO 

^SxP'T 


^Str[a^Ea^S] + ^Str[A<ts + _‘^(^Ui + Dj + Sj)^ - a^V. 


( 11 ) 


The staggered potential, V, splits the tree-level pion masses into hve degenerate gronps, 


= n{mi + rrij) , 


( 12 ) 


according to their S'0(4)-taste irrep, S = I, P,V, A,T.^ It also leads to hairpin (qnark- 
disconnected) propagators with mnltiple poles for flavor-nentral, taste V and A pions. 

We now constrnct the remaining terms in the heavy meson chiral Lagrangian. Ref. [l^ 
showed that, at 0{a^), mixed four-fermion operators with both heavy and light quarks can¬ 
not break taste-symmetry. Because all taste-violation in the Symanzik action comes strictly 
from the light quark sector, discretization errors caused by mixed four-fermion operators 
can be categorized as “heavy-quark errors” and estimated using standard methods QE3- 
Thus the form of the heavy meson portion of the chiral Lagrangian is identical to that in 
the continuum, with the exception that the light quark index can run over both flavor and 
taste. 

Heavy meson chiral perturbation theory (HMyP^ was hrst formulated in Refs. QQ 
and generalized to partially quenched QCD in Ref. 18[- Heavy quark spin symmetry allows 
the pseudoscalar and vector mesons to be combined into a single held which annihilates a 
heavy-light meson: 

1 + ^ 


H. = 




(13) 


where v is the meson’s velocity and a labels the havor and taste of the light quark within 
the meson. Note that, although we use the letter B, the heavy-light meson can be either a 
B, in which the heavy quark is a 6, or a D, in which the heavy quark is a c. We also dehne 
the conjugate field, Ha = which creates a heavy-light meson. 


^ Although we are interested in describing a Euclidean lattice theory, we choose to perform the calculation 
in Minkowski space in order to make intermediate steps comparable to the continuum literature. Our 
results for the form factors will be independent of this choice. 

^ Note that the splitting Ap = 0 because the taste pseudoscalar pion is an exact lattice Goldstone boson 
in the chiral limit. 
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Under heavy-quark spin symmetry, 


(14) 

(15) 


H ^ SH, H ^ HS^ 

S e SU{2), 

while under chiral symmetry, 

H HU\ H ^UH (16) 

U e SU{An). (17) 

Interaction terms between heavy-light and pion helds are constructed using a = = 

exp[i<h/2/], which is invariant under heavy-quark spin symmetry but transforms under chiral 
symmetry as 

a ^ LaU^ = UaR\ = Ua^Ll (18) 

Heavy meson chiral perturbation theory is a joint expansion in the inverse of the heavy 
quark mass, l/mq, and in the residual momentum of the heavy-light meson, k. Thus the 
leading order heavy meson Lagrangian is of 0{k): 

[HaV^ {SabO, + H,] + g^tvn (HaHbY^^Al^) , (19) 

where Vjj = | \a^d^a + = | \a^d^a — ad^a^~\ and tin indicates a trace over Dirac 

spin indices. Combining this with the purely pionic terms, the total chiral Lagrangian for 
heavy-light mesons in which the light quark is staggered is 

/^LO _ /^LO I /^LO 

^ — ^HMxPT “T >t^SxPT* 

III. CHIRAL CORRECTIONS TO B ^ D AND B ^ D* AT ZERO RECOIL 

The hadronic matrix elements for B D^*'^ depend upon six independent form factors: 

{D{v')\cYb\B{v)) = h+{w){v + vY+ h_{v-v'Y, (21) 

{D*{v)\cj^j 5 b\B{v)) = -ihAi{w){w+ l)e*^ 

+ ihA2{u;){v ■ €*)v^ + ihA:i{w){v ■ €*)v'^, 
{D*{v')\c'y>^b\B{v)) = hv{w)e^^''^^e\v'^vg 
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( 22 ) 

(23) 



where w = v ■ v'. In the static heavy quark limit, however, heavy quark spin symmetry 
requires that h_ = = 0 and hj^{w) = = hv{w) = ^(ta), where i{w) is the 

universal function for B decays called the Isgur-Wise function.^ Our goal is to 

calculate the leading nontrivial contributions to the B —> form factors at zero recoil, 

i.e. when v' = v and w = 1. 


At zero recoil, the hadronic matrix elements depend upon only two form factors, h_|_(l) 
and 

{D{v)\cYb\B{v)) = 2n^h+(l), (24) 

{D*{v)\cY75b\Biv)) = -i2e*^hA,il). (25) 


In the static heavy quark limit, ^(1) is normalized to unity corrections to this result 
come from operators of C>(l/mQ). Operators of O^I/tiiq) can be separated into those that 
respect heavy-quark spin symmetry and those that break heavy-quark spin symmetry. The 
former cannot produce logarithmic contributions to 5 —> form factors at leading order 

in yPT because they contribute equally to B{D) and B*{D*) masses, so we do not show 
them here. The single 0{l/mQ) operator that breaks heavy-quark spin symmetry comes 
from the interaction between the chromomagnetic moment of the heavy quark and the light 
degrees of freedom: 


5C = 

mg 


(26) 


This operator generates a splitting between the D and D* meson masses, = {mo* — 
mo) = —A 2 / 8 mc. It also produces a. B — B* mass splitting, but is of 0{l/mb) and can 
be neglected. Finally, we note that there turn out to be no 0{l/mQ) corrections to the form 


factors at zero recoil because of Luke’s theorem [^, so the leading nontrivial contribution 
to h+(l) and hyi^(l) is of 0{l/ml). 

In order to calculate the form factors h+(l) and Ha^ (1), we must hrst map the quark-level 
B —> operator onto an operator in the chiral effective theory: 

C7'‘(l - 75)6 ^ -?(m) Str[foV(l - ( 27 ) 


^ This is true up to radiative corrections, which only affect yPT through a modihcation of the low energy 
constants. 

® Note that the form factor h-{l) appears in the differential decay rate ior B ^ D and is needed in lattice 
determinations of \Vcb\ from B ^ D. Lattice calculations have shown that this term is a small correction, 
and we do not consider here the chiral corrections to this small quantity. 
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B! 

—•- 


\ B 

-*— 


(c) 


FIG. 1: One-loop diagrams that contribute to B ^ D*. The solid line represents a meson contain¬ 
ing a heavy quark, and the dashed line represents light mesons. The small solid circles are strong 
vertices and contribute a factor of . The large solid square is a weak interaction vertex. Diagram 
(a) is a vertex correction; (b) and (c) correspond to wavefunction renormalization. 

We can then calcnlate the desired hadronic matrix elements, Eqs. in the heavy- 

light meson effective theory that is described by the Lagrangian in Eq. CT plus the addi¬ 
tional D — D* mass splitting term, Eq. (ESI). 

The B ^ D and B —> D* matrix elements receive contributions from the diagrams shown 
in Figure d It is necessary, however, to consider these same diagrams at the quark level 
in order to identify sea quark loops. This is because, in SyPT, all sea quark loops must 
be multiplied by 1/4 in order to describe data from staggered lattice simulations in which 
the fourth-root of the quark determinant is taken to reduce the number of tastes per flavor 
from 4 to 1.® Quark flow analysis also allows identification of quark-disconnected hairpin 
diagrams, which can only occur for taste /, V, and A pion loops, that have propagators with 
multiple poles. At the quark level, two vertices appear in the calculation of the B ^ D and 
B ^ D* form factors; they are shown in Figure El The HH'k vertex comes from the LO 
heavy meson chiral Lagrangian, Eq. (P|). and is proportional to The DB vertex comes 
from the weak operator in Eq. EZD- Using these vertices. Figure El shows the same diagrams 
as in Figure m but at the level of quark flow. 

® Throughout this paper we assume the validity of the v^Det trick; for a recent review of the status of the 
v^Det trick see Ref. [^ . 
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at 16 





(a) (b) 

FIG. 2: Relevant vertices at the quark level. Vertex (a) comes from the LO heavy meson chiral 
Lagrangian and is proportional to the coefficient g-,^. The solid line corresponds to a heavy bottom 
or charm quark while the dashed lines correspond to light staggered quarks of any flavor and taste. 
In this vertex, either one or both heavy-light fields must correspond to a vector meson. Vertex (b) 
comes from the B —> D^*') operator in Eq. The solid double line corresponds to the bottom 

quark within the B- or i?*-meson and the solid single line corresponds to the charm quark within 
the D- or Z)*-meson. 



■ - ^ ^ > * 


(c) 

FIG. 3: Quark flow diagrams that contribute to B ^ D and B —> D*. The double line corresponds 
to the bottom quark within the R-meson, the single line corresponds to the charm quark within the 
rW -meson, and the dashed lines correspond to staggered light quarks. Diagram (a) renormalizes 
the R-meson wavefunction while (b) renormalizes the R(*)-meson wavefunction. Diagram (c) 
modifies the R —> D^*'^ vertex. 
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We now present results for the form factors relevant for B ^ D and B D* at zero 
recoil including taste-breaking effects due to the staggered light quarks. 

For the H-l-l-l PQ theory, in which 




1 + 


ml I 16 


Xfi, T.*? ^ \ y\ I / 


j=xu,xd,xs 
H=7,P,4V,4A,6T 




je{MP} 


+ aP' 


V 




dF: 






+ (V^A)i, 


(28) 




1 + 


^a(A) , g. 


+ 


— 1 


1 V—^ 1 

— > F,_ -h - 


ml 487r^/^ 116 

^ j=xu,xd,xs 

E=I,P,AV,AA,QT 




dFxj 

dmxj 




je{MP} 


-I- a 6‘ 


2x' 
V 




- E Df:i({kyZ};{i.v])F, 




+ (f/-kl) , 


(29) 


where x labels the light valence quark within the decaying B^ meson. The first term inside 
the curly braces comes from diagrams with sea quark loops; the flavor index j runs over 
mesons made of one valence quark and one sea quark and the taste index S runs over 
the sixteen pion tastes. The residues and are due to flavor-neutral hairpin 

propagators; their explicit forms, along with the sets of masses are given in the 

Appendix. The second term in the curly braces (with coefficient 1/3) comes from taste- 
singlet hairpins, while the third and fourth terms (with coefficients Pdy and PS'P) come 
from taste-vector and axial-vector hairpins, respectively. The functions F and F are defined 
as 


Fj = F {mj,A^^'>/mj) (30) 

Fj = F {mj,-A^^^/mj) , (31) 
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where 


F {rrij , x) = 


X r A" ^ 3 


Ax + 27r 


\/x^ — l{x‘^ + 2) ^In 1 — 2x{x — V— 1) — ivr'j 


(A'“>)'Mn + 0 |(aW)=>] 


( 32 ) 


and rrij is the tree-level mass of a meson with flavor-taste index j, given in Eq. (IT^ . The 
analytic terms proportional to X+(A) and X^(A) exactly cancel the renormalization scale 
dependence of the F terms. It is interesting to note that heavy-quark symmetry forbids the 
presence of additional analytic terms such as those cx Str(Ad) = (m„ -|- rud -l- mg) or cx a^. 
We have checked that these expressions agree with the continuum partially quenched ones 
when a —> 0 2^. In this limit, the masses of all of the pion tastes become degenerate, so 
Thus the sea quark loop and taste-single hairpin contributions reduce to 
the continuum PQ result, while the taste-vector and axial-vector contributions, which are 
proportional to a^, vanish. 

For the 2-1-1 theory, in which = rUd ^ mg'. 


2+1(1) = 1 + + 


9. 


mt. 




■i — ifni nr-'^i -t'o I- \ yvj / 


J=XU,XU,XS 

E=I ,P, 4 V, 4 A, 6 T 




(5)t. 


Xi 




-|- 5y 








^7 

dFxy, 

dm% 

Xv 


(33) 


r {Bx)PQ',‘2-\-l 


( 1 ) = 1 + 


ml dSvr^/^ I 16 




J=XU^XU,XS 

H=/,P,4V,4A,6T 








<h{<};W})(^)- E BS({<d;W})T2 




+ {V^ A) 


(34) 
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In the case of full (2+1) QCD, the expressions for the residues simplify because QCD is 
a physical (unitary) theory without double poles: 


K 


{Bu)QC D,2+1 = 1 + + 


9^ 


mt 


2 


( 2 -^ 71 = + Fks) rfFni + 


16 


+a^5' 


rriq — ml 

OV TTV 


V 


ml )(m‘i — ) 

nv 11 nv n,,' 


F + 

TTV ' 


2 2 

mt — mi 

Vv ov 


[m. 


w 

2 


I ml, - m2 )(m2 -ml ) 


w 


^Sv ~ 


- hrF —+ (’^ ^ -4) 

[mi —m, m , —mi ) 

^ Vv Vv ^ Vv > / 


(35) 


{S„)QCn,2+1 


. ^ ^+(A) 

(1) = 1 + —iVA + 


mt 


487r2/2 


~ X] (2-^7rs + Fk-s) - T^F^j + -Fr^j 


16 


W5- 


2x' 
V 


^Iv - ^Iv 


9 9 

mz — m% 

7]V w 


_ r j_ 'IV _ r 

(m2^ - m2 ^) (m2 ^ - mj^) (m2^ - i^lv “ ^Iv ) 


w 

ml — mt, \ 

- Iw 2 —rT^+ + (r ^ 71) 

m„ — m , )(m^, — m;! ) / 

1 Vv Vv ^ Vv > / 


TTV } 

^v'v 


Vv' 


(36) 


^(,Bs)QC‘i9,2+l ^ _ 


1 + + 


mt 


167r2/2 


- - A’s, + If,^ 


16 


Ws 


2x/ 
V 


2 2 
mh — ml 

JV ^v 


:Fsv + 


m; 


w 


m. 


TTy 


i^Sv - ^Vv ) i^Sv -^10 - ^Sv ) i'^vv - ) 


?7v 


2 

m , 
’ll/ 


m 


TTV 


(m2, —ml )(m2, — m^ ) 

^ Vv JV'^ riv 


Vv 


+ (V^A) 


Vv 


(37) 




(S. )QCD,2+1 = 1 + + 


9l 


mt 


+a^5' 



487r2/2 
2 


4 E - -Fs, + |f„ 


16 

,2 


me — ml 

oy TTy 


^sv - ^v(^iv 


+ 


9 9 

mz — mz 

rjy TTy 




9 9 

m , — mz 

rjy Fy 


im^, —ml )(m2, —mi ) 

^ Vv ^v'^ Vv 


1-, + (r - 71) 


Tjy 


(38) 


(*) f*) 

Note that there are separate formulae for Bu,d I and Bg ^ Dl in full QCD. This is 


in contrast to the PQ expressions, which are valid for any choice o 
results agree with the continuum full QCD ones when a —0 


light quark flavor. These 
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IV. NUMERICAL ILLUSTRATION OF THE B ^ D* FORM FACTOR 


In this section we present a realistic picture of the behavior of staggered lattice data for 
B ^ D* and compare our S^PT expression to actual staggered lattice data for B ^ D. 

Figure m which shows the full QCD expression for vs. illustrates the im¬ 

portance of accounting for staggered discretization errors in the extrapolation of staggered 
lattice data. There are currently no unquenched staggered lattice data for B D* available, 
so we have added a term linear in to the SyPT expression for and matched onto 

existing quenched data simulated at heavy (> 500 MeV) pion masses P. Thus Figure m 
gives a realistic illustration of what the chiral extrapolation of unquenched data for hA^{l) 
from the MILC coarse lattices (a = 0.125 fm) might look like. The continuum expression 
for has a characteristic cusp at m.,^ = Ac where the internal D goes on-shell. The 

staggered expression has a cusp in the same location due to the taste pseudoscalar pion, 
which receives no taste-breaking shifts to its mass, but the cusp is much milder. 

It is worthwhile to discuss in some detail why the staggered cusp is so mild, or, equiva¬ 
lently, how the staggered curve in Figure |3] becomes the continuum curve when a —0. A 
cusp occurs in Hai every time the internal pion and D go on-shell in the B D* diagram. 
For a staggered pion of taste S, this happens when -|- a^As = A^, where rri^ is the 
tree-level mass of the lattice Goldstone pion and a^A^ is the taste-breaking mass correction. 
Thus, in Figure |3J there is a cusp in the staggered curve every time = A^ — a^As. On 
the MILC coarse lattices, all of the 0(a'^) mass-splittings are greater than Ac = 0.14 GeV, 
so the additional heavy staggered pions do not produce cusps in Figure El The single stag¬ 
gered cusp due to the lattice Goldstone pion is small because it is weighted by 1/16 (from 
the average over pion tastes in the loop) as compared to the continuum one. As the lattice 
spacing is reduced, more and more tastes will be able to produce cusps to the left of the 
continuum one. These cusps will begin at = 0 and move to the right as the lattice spac¬ 
ing becomes smaller. Finally, at a = 0, all of the cusps from the non-Goldstone pions will 
come to rest at the location of the continuum one, and the sum of these sixteen staggered 
cusps will equal the single cusp in the continuum curve. In addition to softening the cusp, 
the heavy staggered pions decrease the curvature due to chiral logarithms in this is a 
generic effect of taste-breaking. Thus the staggered data are expected to be almost linear, 
even when the continuum result is not. 
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0.925 



0.895 


0 


0.1 

m/ (Gev') 


0.2 


FIG. 4; Qualitative behavior of h^{l) vs. m^. The overall linear contribution comes from 
matching to existing quenched data |^. The curve with the large cusp is the continuum expression, 
whereas the (dashed) curve with the mild cusp includes staggered discretization effects. We use 
the measured values of the pion mass-splittings and taste-breaking hairpins from the MILC coarse 


lattices as input into the staggered curve 


s 


In practice, one extrapolates staggered lattice data to the continuum by first fitting 
to Eq. dHEI) and then removing taste-breaking discretization errors by setting the terms 
proportional to a? in Eq. dSED to zero. We note that simulations are not likely to be sensitive 
to the cusp anytime soon, even if staggering did not smooth it out, because the cusp only 
occurs at values very close to the physical pion mass. Thus, in the case oi B ^ D* ^ it is 
especially important to use S^PT to extrapolate to the physical light quark masses. 

We can also directly apply our SyPT expression to the available unquenched data for 
5 — > H j^. Figure m shows h+(l) vs. in full QCD. In this case, the pion and D* in 
the loops of the diagrams of Figure Q cannot go on-shell, so there is no cusp as in the case 
oi B ^ D*. The dashed line is the result of a ht of the staggered expression to the three 
data points. The solid line is the continuum extrapolated curve, while the square is the 
continuum extrapolated value of h+(l) at the physical value of the pion mass, with error 
bars. The difference between these curves is very small, and the extrapolated value hardly 
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FIG. 5: /i+(l) vs. m\. The three full QCD data points (circles) were calculated on the MILC 
coarse lattices (o = 0.125 fm) j^. The upper (dashed) curve is a fit to the data using the complete 
staggered formula, while the lower (solid) curve is the continuum extrapolated curve. The square 
is the extrapolated value of /i+(l) at the physical pion mass with error bars. 

differs from the result of a naive linear fit. Nevertheless, the S^PT analysis is useful for 
this quantity because it demonstrates that the systematic errors associated with the chiral 
extrapolation are small. 


V. FINITE VOLUME EFFECTS IN B ^ UP) 

The functions F and F, which appear in h+(l) and h^^(l), respectively, are modihed 
by the hnite spatial extent of the lattice. Using the formulae for hnite volume corrections 
to typical HMyPT integrals given in Ref. we hnd that F{m,A) receives the following 
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correction due to the finite lattice volume: 


(5Fpv(m, A,L) = 


E( 


m 


n^O 


\128xy^ 


e-y) 


4™10 


y^x — {y — lQ)y X 


+2{y + 2A)y'^x^ — 16{y — 2)y^x^ + QQy'^x'^ — 128y’^x^ + 256?/ 


-a/^ 


y'^x^^ + {y — 15)?/^x^ — (3?/ + 35)y‘^x'^ 


+ {16y^ - 27y + 15)yx^ + i-112y^ + lly - 9)x^ + 256y^x 

y^x^^ — (y — 16)y^x^^ 

+2(y + 24)yx^ — 16(y — 2)yx^ + ^Qyx'^ — 128yx^ + 256 


—2567ry^ — vry^e'^^^'^^erf ^ ^ 


(39) 


where x = A/m as before, y = nmL, and n = \/W. The correction to F is identical except 
for X — > —X. This formula was derived as a series expansion in l/(nmL). In our numerical 
evaluation of this formula we truncate the sum to the values of n = 1, a/2, a/3, a/4, a/5 and 
a/ 6T An expansion in x = A/m shows that the leading contribution to 5Fpy is proportional 
to A^, as expected: 

^F^y(m,A,L) = ^ 

Figure ini shows the contribution to h+(l) in full QCD from finite volume effects for the 
MILC coarse lattice (a = 0.125 fm, L = 2.5 fm). Recall from Figure El that h+(l) is close 
to one, whereas the finite volume corrections in Figure El are less than 10“*^ in the range 
of pion masses relevant for current staggered lattice simulations. The size of finite volume 
corrections to h^/l) are similarly small. We therefore conclude that finite volume errors 
are negligible in both the B ^ D and B D* form factors, and can be accounted for as 
an overall systematic error in lattice calculations, rather than subtracted before the chiral 
extrapolation. 
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determined that truncating the sum at n 


mL > 2.5. 


•a/ 5 approximates the full answer well (~ 3%) for 
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FIG. 6: Finite volume correction to /i+(l) as a function of Recall that /i_|_(l) is close to 1, so 
these corrections are smaller than one part in 10“'^ in current staggered simulations. 

VI. CONCLUSIONS 

In this work we have calculated the B ^ D and B — D* form factors at zero recoil 
to NLO in SxPT. We have presented expressions for both a “1+1+1” partially quenched 
theory {niu ^ rrid ^ mg) and a “ 2 + 1 ” partially quenched theory (m„ = rrid 7 ^ m^,), as well 
as for full (2+1) QCD. These formulae apply to simulations in which only the light quark 
is staggered. They include 0(0^) taste-breaking discretization errors, and are necessary 
for correct continuum and chiral extrapolation of staggered B D^*'> lattice data. Use of 
these expressions, along with the double ratio method of Ref. P] and in combination with 
experimental input, should allow a precise determination of the CKM matrix element \Vcb\- 
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APPENDIX 


In this section we collect formulae necessary for understanding our form factor results. 
We follow the notation of Ref. 

The residues and appear because of single and double poles, respectively, in 

the flavor-neutral hairpin propagators: 


.,/c] 


d 


2 ^’ 


{a}) ^ {m}). 


(Al) 


Once one takes the mass of the overall flavor-taste singlet pion (which corresponds to the 
physical r]') to infinity, the relationships among the taste-singlet pion masses simplify: 


= '"t, = ™D,. 
2 , 2m|, 


(A2) 


Thus the following mass combinations appear in the H-l-l-l (m„ 7 ^ rud 7 ^ PQ result: 

{/i} = {mu,mD,ms} (A3) 

When the up and down quark masses are degenerate, the pion mass eigenstates become: 

,2 „„2 „„2 

V/’ 

3 2 e/ 

-a 6y 



= mu, 

V=^Dy, 



mly + ml^ 



^Uv + m% 


3 

4^ 




Z = 




a25: 


V 


(^Sv-Kv) + 


9(a2h(,)2 


2 ^ Uv' X 6 

Thus the following mass combinations appear in the 2 - 1-1 (m^ = 7 ^ m*) PQ result: 


(A4) 


{M^x) = {m^.rnx}, 

{M^x) = {m^,rn^i,mx}, 

{p} = {mu,ms] (AS) 
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